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SUMMARY 

A simplified stochastic model is proposed for crack initiation and 
short-crack growth under creep and creep-fatigue conditions. Material Inhomo- 
geneity provides the random nature of crack initiation and early growth. In 
the model, the influence of microstructure is introduced by the variability of 
(1) damage accumulation along grain boundaries, (2) critical damage required 
for crack initiation or growth, and (3) the grain-boundary length. The proba- 
bilities of crack initiation and growth are derived by using convolution inte- 
grals. The model is calibrated and used to predict the crack density and 
crack-growth rate of short cracks of 304 stainless steel under creep and creep- 
fatigue conditions. The mean-crack initiation lives are predicted to be within 
an average deviation of about 10 percent from the experimental results. The 
predicted cumulative distributions of crack-growth rate follow the experimental 
data closely. The applicability of the simplified stochastic model is dis- 
cussed and the future research direction is outlined. 


INTRODUCTION 

Life prediction is an important parameter in evaluating the safety and 
reliability of structural components for high-temperature applications such as 
rocket engines, gas turbines, and nuclear powerplants. Crack initiation and 
early growth constitute most of the life of components, especially under creep 
and creep-fatigue conditions. Experimental observations (Ohtani et al . , 1983, 
1986, 1987) have shown large fluctuations in the crack initiations and growth 
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rates of short cracks. These fluctuations are caused by th a randomness in the 
grain size, the local stress, and the resistance to local failure. Few analy- 
ses have been conducted of the stochastic nature of the growth of microstruc- 
turally short cracks. However, a considerably larger number of investigations 
(Kozin and Bogdanoff, 1981; Virkler et al . , 1979; Lin and Yang, 1983; Ishikawa 
et al . , 1987; Spencer and Tang, 1988) have been carried out for long-crack 
problems. Long cracks usually show smaller fluctuations in the crack-growth 
rates than short cracks. The two approaches taken for short cracks thus far 
are based either on Monte Carlo simulations (Kitamura and Ohtani , 1987, 1988) 
or on the randomization of empirical crack-growth equations (Cox and Morris, 
1987). The approach taken here is based on the damage accumulation along grain 
boundaries and the critical damage required for failure under creep and creep- 
fatigue conditions. 1 The fluctuation of the crack-growth rate is formulated 
with the fundamental reliability analysis (Thoft-Chri stensen and Baker, 1982), 
which is based on the variability of the damage accumulation and the critical 
damage of the grain boundaries. The model is calibrated and is then used to 
predict the crack initiation and early growth rate of short cracks under creep 
and creep-fatigue conditions. 


STOCHASTIC MODEL FOR CRACK INITIATION AND GROWTH 

For many polycrystalline alloys under creep or creep-fatigue conditions, 
cracks initiate and grow along preferential grain boundaries. This localiza- 
tion of the cracking sites is due to the operative high-temperature failure 
mechanisms. Growth and coalescence of cavities and grain-boundary sliding are 
examples of local failure mechanisms leading to grain-boundary cracking (Garo- 
falo, 1965). The time required for the initiation of a crack differs between 
seemingly identical grain boundaries because of creep resistance and geometric 
differences. The crack-growth rate is also dependent on such things as the 
grain-boundary orientation with respect to the loading axis, the distance 
between triple points, and the distribution of the precipitates and ledges on 
the boundaries. The highly random crack initiation and growth processes are 
created by the presence of mi crostructural inhomogeneities. The local loading 
conditions along grain boundaries are also random because of creep deformation 
anisotropy, grain shape, and constraints. A simplified stochastic model is 
proposed for predicting the probability of crack initiation and growth under 
creep and creep-fatigue conditions. This model takes into consideration mate- 
rial inhomogeneities, microstructure geometric variations, and local loading 
randomness . 

The assumptions included in the model are that (1) cracks grow only along 
grain boundaries, (2) crack length is measured as the projected length on a 
plane perpendicular to the loading axis, (3) projected grain-boundary length 
between adjacent triple points d is a random variable with density function 
f(d), and (4), the width of the crack is ignored. The two-dimensional model 
can be modified, if necessary, to a three-dimensional model by converting the 


Hhe term "creep fatigue", as used here, implies a reversed cyclic creep 
deformation at high temperature under which creep dominates the failure proc- 
ess. A discussion of synergistic creep-fatigue interaction is beyond the scope 
of this report. 


2 


crack length into crack area. It is also assumed that crack initiation and 
growth occur in a discrete manner, increasing in segments that are equal to 
the distance d between two triple points after a certain time t (fig. 1). 

It has been observed experimentally that grain-boundary triple points act as 
crack arresters (Ohtani et al . , 1984). Therefore, crack growth is usually 
halted at triple points or grain-boundary kinks (sharp bends) until enough 
damage is accumulated in the adjacent grain boundary for the crack to extend 
again. It is postulated that crack initiation and growth occur when the accu- 
mulated damage <t> equals a critical value 4> c . The critical damage <J> C can 
be interpreted as the resistance of a particular grain boundary to cracking, 
which is a random variable. The density function of the critical damage $ c 
is given by g(<j> c ). The damage accumulation <t> is the local driving failure 
parameter for a particular grain boundary. The critical damage <t> c is inde- 
pendent of the applied stress but 4> is strongly related to it. The damage 
accumulation <j> is the integral of the damage accumulation rate <j>: 


nt 


4, = 


<j> dt 

J 0 


( 1 ) 


The damage rate <j> is a random variable that is a function of the local 
stress ai . The local stress state ai is a function of the constant global 
applied stress a g and the crack length 2. Assuming a steady-state failure 
process, the damage accumulation rate is taken to be independent of time. 

♦ - <j>t (2) 

The density function k of <i> is related to the density function of <j>, 
h(<j> | ag , 2.) by 

k(4>|o g ,2) = 1/t • h(<j>|a g ,2) (3) 

If the grain boundary under consideration is adjacent to a crack, the failure 
is called crack growth. If the grain boundary is isolated from other cracked 
regions, the failure is called crack initiation. 

This model is applicable, but not limited, to monotonic creep conditions. 
It can also be applied to time-dependent fatigue (i.e., creep-dominant 
fatigue). As an example, under a slow-fast fatigue loading (fig. 2) of 304 
stainless steel at 973 K in vacuum, irreversible grain-boundary sliding was 
observed. This was due to the creep deformation in the tensile half of the 
cycle (Taira et al . , 1978). Under those conditions, the intergranular crack 
initiation is associated with creep damage accumulation in the tensile half of 
the cycle. Thus, the time should be converted into the number of cycles and 
the local stress ai should be interpreted as an equivalent stress of the ten- 
si.le half of the loading cycle. 

This damage accumulation formulation is expressed in mathematical terms 
only. Physical explanation of the damage equation is discussed next. If one 
considers that the damage is associated with cavity growth, the damage function 
4> can be written in terms of the number and the radii of the cavities 
observed on a particular grain boundary. The damage function <J> can be formu- 
lated as the ratio of the summation of cavity length over the grain-boundary 
length. Experimental observations have shown that the summation of the length 
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of grain-boundary cavities is proportional to time (Yang, 1984). Hence, a 
physical meaning is given to the steady-state damage accumulation of 
equation (2). More sophisticated damage equations can, of course, replace 
equation ( 2 ), which can be written in terms of cavity initiation and growth 
1 aws . 


PROBABILITY OF CRACK INITIATION AND GROWTH 
Crack Initiation 

As seen previously, failure on a specific grain boundary occurs when the 
damage accumulation exceeds a critical damage level. Let us now define a new 
variable <|> m , known as the failure function, as the difference between 4 > c 
and <j>: 


♦m a ♦c ~ ♦ (4) 

The failure criterion is simply reduced to <j> m = 0 . Consequently, the proba- 
bility of a crack initiating at a particular grain boundary Pf(t) is the 

probability that <j> m is less than or equal to zero: 

Pf(t> = P< 4 > m < 0 ) (5) 

The density function of the new random variable <j> m is given by m 0 (<j>Jt, a g ) 

for a given time t and global stress a g . The density function m 0 is 

related to the density functions of the two independent random variables <|> c 
and <j>, by the following convolution integral: 


m (<j>| t,o„> 
o T m 1 g 


oo 

9 % + 4 >> k ( 4 > [ t , cr g , a = 0 ) d<J> 

^0 


Knowing that 41 = <j>t, then 


( 6 ) 


m (<|> m | 
0 T m' 




<3% + *t> h<«j»|a , i = 0 ) d<j> 

J 0 


(7) 


If the density functions of 4 > c and <j> are known, then the probability Pf(t) 
for crack initiation at one particular grain boundary before time t for a 
given global stress o g is given by 


p.(t) = P(<b < 0) 
f T m - 


ITL<$J t,a_) d«j)_ 
0 m 1 g m 


( 8 ) 


A polycrystalline material has many grain boundaries, which are potential 
crack initiation sites. The probability that n number of cracks initiated 
in time t for a total number of potential crack initiation sites N is 
given by the following binomial distribution: 
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( 9 ) 


N! [P,(t)] n [1 - P f (t)] N ' n 

(P ) _ 1 1 

' f n (N - n)! n! 

The expected value or mean number of crack initiations in a given time t and 
its variance are given respectively by 

E(n) = N P f (t) (10) 


and 


VAR = E([n - E(n)] 2 ) = N P f (t) [1 - P f (t>] 


( 11 ) 


Crack Growth 


The criterion for crack growth in terms of accumulated damage is similar 
to that for crack initiation. The crack-growth criterion still corresponds to 
the damage accumulation <J> reaching a critical value, <J> C ; however, the damage 
accumulation rate is accelerated by the presence of a crack. The failure func- 
tion <|> m of equation (4) for the grain boundary i immediately adjacent to 
the crack tip is redefined for crack growth as 


4»m = 4>c “ $i(fl-»og)t (12) 

where 4> i ( 2. , ag ) is the damage accumulation rate at grain boundary i, and t 
is the time elapsed after the crack reaches the particular grain boundary i. 
Since crack growth occurs when <j> m = 0, the time interval ti required for the 
crack to travel the entire grain boundary i is given by 

ti = 4 > c/ ( j > i (13) 


The density function of the time t, in terms of the density functions g(<j> c ) 
and h(<j>j) can be shown to be 


s(t. |o f ,a> 


<i). g(<t> i t- | Og) h(<j>. |a g ,&) d<}>. 

^0 


(14) 


The crack-growth rate of a short crack is then given by 

da/dt = di /ti (15) 

where di is the projected distance between two triple points of grain bound- 
ary i. The density function of the crack-growth rate in terms of the density 
functions s ( t i ) and f ( d i ) of t, and d^ , respectively, is given by 


co 

P 


P (dft/dt|a ,ft) = 



f (da/dt 


V S(t il a g> a) dt i 


(16) 
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The damage accumulated before the grain boundary i becomes adjacent to the 
crack tip is assumed negligible because it was shown experimentally that pre 
damage of stainless steel specimens has little effect on crack growth under 
creep and creep-fatigue conditions (Ohtani and Kitamura, 1986). A formulation 
that takes into consideration the predamage effect is described in the 
appendix. 


Normal ization 

As seen previously, the damage accumulation rate <j> is a function of the 
local stress a], which depends on the applied global stress a g and crack 
length Q.. For simplicity, the following relation can be assumed: 

a] = K(u|!l)ag (17) 

where K(u|&) is a stress concentration factor function of the crack length 2. 
and mi crostructural randomness u. Under creep conditions, the damage accumu- 
lation rate can be taken to be a power function of the local stress 

<j> = Ao“ (18) 

where A and a are assumed to be material constants. Assuming a steady- 
state damage rate, the damage accumulation is then given as a function of the 
global stress and K( u | fi.) by substituting equations (17) and (18) in 
equation (2): 


<(> = <j>t = A(K(u|a)a g ) a t 


(19) 


where 4> = A K(u|2) a 


4> = 4>t 



(20) 

( 21 ) 


Here, <j> is a random variable because K(u|2) is a random variable, but t is 
a deterministic variable. The density function fi of the normalized damage 
accumulation rate is related to h of <j> by 


fi(<j>) = o“ h(<j>) 


( 22 ) 


The probability of the crack initiation can be rewritten in terms of the nor- 
malized variables as 


P f (t) 


r ° 

fiL(<b|f) d* = o®P- 
o m 1 T m g f 

J -00 


(23) 


where 


m Q (<l>Jf) = 9(<t> m + R(<H& = 0) d<j> 


(24) 
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The density function of the crack-growth rate in terms of the normalized time 
interval l\ = ogt} can be shown to be equal to 


P (da/dt|&> = 


r°° 

l. f (dft/df 

^0 


t. ) s(^ |a) df . 


(25) 


scf, | a> = o® s(t i |o ,a> (26) 

Therefore, the actual density function is related to the normalized density 
function by 

P (da/dt|<T ,&) = (l/cr ) a P „(da/df|&) (27) 

g ' g g g 


It should be noted that Pf and Pg are independent of a g . The terms Pf 
and Pg are easily calculated for any arbitrarily applied global stress og 
by equations (23) and (27), respectively. 


APPLICABILITY OF STOCHASTIC MODEL 

Few experiments have been carried out to determine crack initiation and 
growth of short cracks at high temperature. Because of the small data base 
available, only the applicability of the stochastic model Is discussed. This 
model is applied to two different data sets for which experimental data are 
available. Both experiments were carried out on smooth specimens of 304 stain- 
less steel with an average grain diameter of 40 pm. The cracks initiated and 
grew on the specimen surface. In the first experiment (Ohtani et al . , 1983), 
the number of crack initiations was monitored under monotonic creep at 923 K 
in air at applied stress levels of 98.1 MPa and 147.1 MPa. In the second 
experiment, crack initiation and growth were measured under slow-fast fatigue 
(Ohtani et al . , 1986). The total strain range was equal to 1 percent and the 
strain rates in tension and compression were 10 -3 percent/s and 1 percent/s, 
respectively, at 923 K in vacuum. 

For the model calibration, the density functions of the damage rate 
h(<t>|ag,Sl) and of the critical damage g(<t> c ) should be determined from actual 
experimental data. However, standard density functions are assumed in this 
study because of the limited data available. Since 4> c and <j> are positive 
variables, the two-parameter Weibull and logarithmic normal distributions are 
adopted here for 4> c and <t>, respectively, and are given below for 
completeness: 


9 < 4 > c > 




(28) 
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h(<j> | a , Q.) = e (29) 

j V s ’in** 


where |3 and b are material constants, nty is the mean value of <j>, 
and a] n ^) is the standard deviation of the In <j>. The normalized function 
of the damage accumulation rate Fi($|!l) and the density function of the damage 
accumulation k($jcg,t,&), which are derived from equation (29), are given by 


fi(<fc|H> 


yjlv 


— e 


1 n<}» q 


In [*- 


ln<j> 


k<4»| <y n . t ,ft> = 


1 


mjt 

v i 


ln<j> 


"In** 


( 30 ) 


( 31 ) 


where m^ is the mean value of <j>. 

The exponent a, in equation (23), is given by (y + 1) where y is the 
creep exponent of the material. The choice of a to be equal to (y + 1) is 
based on the assumption that the damage accumulation is proportional to the 
local strain energy density rate U®eoi , which is proportional to 
for a power-law creep material. The creep exponent y of 304 stainless steel 
is equal to 7 under monotonic creep conditions at 923 K in air (Kitamura and 
Ohtani , 1988). 


Crack Initiation 

The constants in the density functions of the critical damage and the 
damage rate are calibrated with actual experimental data under monotonic creep 
at 98.1 MPa. The constants are adjusted for the theoretical expected value of 
the number of crack initiation sites E(n) to fit the actual data (Ohtani 
et al . , 1983). The total number of potential crack initiation sites N is 
approximately equal to 900 sites/mm 2 for this material. The constants obtained 
from the calibration are listed in table I. The expected value of the crack 
initiation sites is plotted versus time and compared with the experimental 
results (fig. 3(a)). The rate of change of the expected value E(n) increases 
with increasing time. Also shown is the scatter band of the mean value, plus 
and minus one standard deviation. The calibration fits the data well in view 
of the restriction of the assumed density functions for 4> and <t> c . 
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As a check to the validity of the calibration, the expected value of the 
number of cracks inti a ted is calculated for a different stress level and then 
compared with experimental data. Figure 3(b) shows the predicted results under 
creep conditions at an applied stress level of 147.1 MPa. The experimental 
results are shown to be within one standard deviation and are always lower than 
the predicted mean value E(n). The stochastic model predicts the mean-crack 
initiation lives to be within an average deviation of about 10 percent. 


Short-Crack Growth 

The stochastic model for early crack growth is applied to 304 stainless 
steel under creep-fatigue condition at 923 K (Ohtani et al., 1986). The con- 
stants fi and b of g(4> c ) are assumed to have the same values as in the pre- 
vious crack initiation study under creep condition because the variability of 
the critical damage is only material dependent. The acceleration of the damage 
accumulation rate h(<j>), caused by the presence of a crack, is taken to be a 
simple linear variation of the mean with crack length 2 as given by 

m<j> - Co + Ci 2 (32) 

The standard deviation is assumed to be independent of 2 for simplicity. 

The constant Cq corresponds to the mean value of <(> for 2 = 0, which can 
be determined easily from crack initiation data as described in the preceding 
section. The constant C] is calibrated by using long crack-growth rate data 
of notched specimens with crack length greater than 1 mm. Note that Ci is 
calibrated by using long-crack notched specimens that are different from the 
smooth, short-crack specimens. The calibrated constants Cq and C] are 
listed in table II. The distribution of the grain-boundary length f(d) is 
calculated by assuming a normal distribution with an average grain length of 
0.02 mm and a standard deviation of 0.005 mm. 

The sensitivity of the damage accumulation prior to the dominant crack's 
reaching a particular grain boundary is investigated first with the formulation 
given in the appendix. The values of the constants C] and Cq used in this 
study are listed in table II. The analysis reveals the probability that the 
crack-growth rate d2/dt is insensitive to the predamage time t Q for 
2 2 0.03 mm. When the crack reaches the grain boundary under consideration, 
the damage accumulation due to stress concentration is accelerated to a much 
faster rate than the predamage rate. This is due to the large difference 
between the values of C] and Cq. These calculations are consistent with 
experimental observations (Ohtani and Kitamura, 1986) indicating that, once a 
crack appears, the overwhelming damage is localized in the grain boundary imme- 
diately adjacent to the crack tip. The predicted cumulative probability dis- 
tributions of short crack-growth rates for two different crack lengths of 
0.03 mm and 0.06 mm are given in figure 4. The experimental results are shown 
for purposes of comparison. For the given density distributions of <j> c and 
<j>, the stochastic model gives the same range of crack-growth rates as the 
experimental data. The predicted distributions show good correlation with the 
experimental results for higher crack-growth rates. For lower crack-growth 
rates, the small disagreement may be due in part to the difficulty of monitor- 
ing extremely slow crack-growth rates. The dependence of the crack-growth rate 
on crack lengths is plotted in figure 5. The mean crack-growth rates with the 
10 percent and 90 percent confidence lines are shown for crack lengths ranging 
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from 0 . 03 mm to 0.2 mm. Also shown are the experimental upper and lower crack- 
growth rates cf appi oximately 50 initially monitored cracks. The predicted 
mean value falls between the upper and lower values of the experimental crack- 
growth rates. The predicted 90 percent and 10 percent confidence lines follow 
closely the experimental upper and lower limits, respectively. 


SUMMARY OF RESULTS 

A simplified stochastic model was proposed to predict the distribution of 
the initiation of cracks and their early growth under creep and creep-fatigue 
conditions. The model was formulated by using concepts of damage accumulation 
and critical damage required for grain-boundary failure. The density functions 
of the damage rate and the critical damage were assumed to have standard forms 
(log-normal and Weibull distributions, respectively). The mean and the vari- 
ance were calibrated by comparing the predicted results with experimental data. 
Damage history and the presence of a dominant crack were introduced in a sim- 
ple form that influenced the mean of the damage accumulation rate. The stan- 
dard deviation of the damage rate was assumed constant in this study. The 
model predicted mean-crack initiation lives at different stress levels within 
an average deviation of about 10 percent. The cumulative distributions of the 
crack-growth rates under creep-fatigue conditions for different crack lengths 
followed closely the variability of the experimental data. 

Although the proposed model gave good predictions, there is a need for 
model improvements and sensitivity studies of the damage distributions and 
their dependence on crack length and predamage. The damage accumulation should 
be formulated from a micromechanics, or cavity growth, law. The crack-growth 
behavior was oversimplified by allowing for only one dominant propagating 
crack. Multiple cracks and crack coalescence, based on the formulation in the 
appendix, should be considered. Damage recovery must be accounted for in the 
formulation for more complex loading conditions with compressive segments. The 
effect of predamage formulation in the appendix must be verified with additional 
experiments . 
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APPENDIX - PREDAMAGE EFFECT ON CRACK GROWTH 


The failure function 4> m of equation (4) with predamage is now expressed 
as 


= (j> c - <j> 0 i - <j>i(&.tfg)t (Al) 

where ^(a.ag) is the damage accumulation rate at boundary 1 when the crack 
is immediately adjacent to that particular boundary, t is the time elapsed 
after the crack reaches the particular grain boundary i, and 4> c is the crit- 
ical damage needed for the grain-boundary failure. The function <b 0 \ is the 
damage accumulated at that grain boundary prior to becoming adjacent to the 
crack tip. In other words, <{> 0 j represents the damage state of the grain 
boundary prior to its being located adjacent to the main crack. Assuming that 
the dominant crack affects only the adjacent grain boundary, the damage 4> 0 i 
can be expressed as a function of the damage accumulation rate: 


i-1 


*oi ■ * 0) ‘j 


- Vv 1 - 0) *o 


(A2) 


where tj is the time interval for the crack to travel the entire length of 
grain boundary j. Summation of tj is equal to the total time t 0 elapsed 
before the crack reaches the grain boundary i. The density q(<j>^) of 
(<|>oi + 4>i t > is given by: 


q(<l> t |o g ,a,t,t 0 > = 


oo 

p 




(A3) 


where h 0 i ($} (og.ft = 0), <j>i (ag.Q,)) is the joint probability density of 
<j> j ( on , 2. = 0) and <j>i(ag,2.). The density function m(<j> m ) of the failure 
function is therefore given by 


m < <t> m kg,a,t, t o ) 


J 0 


g(<t> m + 4> t > q<* t |a g ,H,t,t 0 > d* t 


(A4) 


The cumulative probability function S< t i |ag,a,t 0 ) of the time interval 
required for the crack to travel the entire length of grain boundary i is 
given by 


nO 


set. 


IV a ’V s 


m(< ^m’ a * t i* t o ) 


d* 


m 


(A5) 


The density function of the crack-growth rate dd/dt is given by equation (16) 
with s ( t i |ag,Q.,t 0 ) being the derivative of equation (A5). However, a special 
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constraint must be placed on the function s < t i ) . It is assumed that no crack 
growth occurs on a particular grain boundary unless that boundary is immedi- 
ately adjacent to the main crack. Therefore, s ( t i ) for t^ < 0 must be zero. 

The following two requirements are needed to satisfy the positive constraint 
on s ( t i ) for ti =0: 

(1) The function <|> 0 i does not exceed 4> c at grain boundary i before 
the crack reaches that grain boundary. 

(2) If <j> 0 i is greater than <j> c at grain boundary i, the boundary 

acts as a potential crack. The main crack instantaneously jumps to the next 
boundary, i + 1, when it reaches grain boundary i. 

Case (1) is consistent with the assumption shown in figure 6(a) for a sin- 

gle crack growth. Case (2) gives future directions for handling crack coales- 
cence for multiple crack growth (fig. 6(b))'. The procedure for adjusting 
s(tj) to be always positive for case (1) is similar to the adjustment of the 
density function of strength after proof-testing (Brent-Hall, 1988). The 
adjusted function, shown schematically in figure 7(a), is given by 



For case (2), the negative part of ti is lumped in a delta function at 
0 + with a magnitude equal to S( t i ) . The adjusted density function is then 
given by r 


s'^log.a.y - 


t. < o 


s ( t . ) + s«n s(t 1 =o) t 1 > o 


This relation is illustrated in figure 7(b). The density functions Pg, for 
cases (1) and (2), are schematically drawn in figure 8. For both cases, the 
crack growth accelerates as the historical time t 0 increases. The infinite 
crack-growth rate, dd/dt = ® (fig. 8(b)), corresponds to the coalescence of the 
main crack with a potential crack. Crack coalescence is an important phenome- 
non that can shorten creep life drastically. However, the exact solution is 
beyond the scope of this report because a multiple crack-growth model is 
required. 
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TABLE I. - CALIBRATED CONSTANTS FOR DISTRIBUTIONS OF <t> c , AND <(> 
FOR CRACK INITIATION ANALYSIS FOR M0N0T0NIC CREEP OF 
304 STAINLESS STEEL 


Constants for 
density function 
of critical 
damage (Eq. 28), 

Constants for density 
function of damage 
accumul ated ( Eq . 29) , 

m) 

Constants for normalized 
densi ty function of 
damage accumulated 
(Eq. 30), 

E(4>) 

B = 2 
b □ 10 

°1 n<|> ° ’-0 4 

nv = 3.4x10 
(o g =98.1 MPa) 

°ln4> = 1>0 
" 3.9xl0- 20 


TABLE II. - CALIBRATED CONSTANTS FOR 
DISTRIBUTIONS OF 4> c AND 4> FOR 
CRACK-GROWTH ANALYSIS UNDER 
CREEP-FATIGUE CONDITIONS OF 
304 STAINLESS STEEL 


Constants for 
density function 
of critical 
damage (Eq. 28), 

g(4- c > 

Constants for density 
function of damage 
accumulated for crack 
growth as a function 
of crack length a 
(Eq. 29), 

h(4>|og,6) 

0 = 2 

m l n<J> = C 0 + C l fi 

b = 10 

C 0 = 0.008 


C, = 1.08 


°1 n4> = l-° 
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(a) PROJECTED GRAIN-BOUNDARY LENGTH, dt 








FIGURE 2. - STRAIN VARIATION FOR SLOW-FAST CREEP FATIGUE OF 30M 
STAINLESS STEEL. £ c = STRAIN RATE IN COMPRESSION; £ t = STRAIN 
RATE IN TENSION; £ t « E c . 


(b) CRACK INITIATION. 



(C) CRACK GROWTH. 

FIGURE 1. - SCHEMATIC VIEWS OF CRACK INITIATION 
AND GROWTH. 
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CUMULATIVE PROBABILITY 




TIME, t, h 

(a) CALIBRATED. GLOBAL STRESS LEVEL, 98.1 MPa. (b) PREDICTED. GLOBAL STRESS LEVEL, 147.1 MPa. 

FIGURE 3. - CRACK-INITIATED DENSITY IN MONOTONIC CREEP OF 304 STAINLESS STEEL. 923 K IN AIR. E(n) = MEAN NUMBER OF INITIATED 
CRACKS; 0 = STANDARD DEVIATION 




CRACK-GROWTH RATE, di/dt , mm/CYCLE 

(a) CRACK LENGTH, i = 0.03 mm. (b) CRACK LENGTH, / = 0.06 mm. 

FIGURE 4. - CUMULATIVE DENSITY FUNCTIONS OF CRACK-GROWTH RATE OF 304 STAINLESS STEEL FOR TWO CRACK LENGTHS UNDER CREEP- 
FATIGUE CONDITION. 923 K IN VACUUM; SLOW-FAST FATIGUE; STRAIN RANGE, A£, 1 PERCENT. 
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FIGURE 5. - COMPARISON OF PREDICTED AND EXPERIMENTAL VARI- 
ATION OF CRACK -GROWTH RATE OF 30tJ STAINLESS STEEL AS 
FUNCTION OF CRACK LENGTH UNDER CREEP-FATIGUE CONDITION. 
923 K IN VACUUM; SLOW-FAST FATIGUE; STRAIN RANGE, A£, 

1 PERCENT. 


t = t. 


GRAIN BOUNDARY i 
DAMAGE ACCUMULATION 



CRACK GRAIN BOUNDARY j 

INITIATION DAMAGE ACCUMULATION 


M 

CRACK GROWTH 


t = to 


CRACK GROWTH 


CRACK INITIATION 
AND GROWTH 



CRACK GROWTH 



CRACK COALESCENCE 


(a) FOR SINGLE CRACK GROWTH. 


(b) FOR CRACK COALESCENCE. 


FIGURE 6. - SCHEMATIC VIEWS OF CRACK-GROWTH PROCESS. t 1 < t 2 <t 3> 
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DENSITY FUNCTION, Pg 



(a) FOR SINGLE CRACK GROWTH. <b> FOR CRACK COALESCENCE. 

FIGURE 7. - DIAGRAMS OF S (tj ) AND s' (tj ) DENSITY FUNCTIONS OF TIME INTERVAL tj. THE IMPOSSIBILITY OF CRACK 
GROWTH BEFORE tj = 0 IS REPRESENTED BY CORRECTED DENSITY FUNCTION S'. 



(a) FOR SINGLE CRACK GROWTH. (b) FOR CRACK COALESCENCE. 

FIGURE 8. - DIAGRAM OF DENSITY FUNCTIONS OF CRACK-GROWTH RATE. GLOBAL STRESS, Og, CONSTANT; CRACK LENGTH, l. 
CONSTANT; PREDAMAGE TIME, t Q ; t Ql < \ q2 < t o3 . 


l 
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